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Chapter I 
INTRODUCTION 

Intenlaminar -stresses play an Important roTe in the load trans- ’ 
fer mechanism In composite laminates. Both numerical and experimental- 
results have demonstrated that when a thin laminate is subjected to a 
uniaxial extension (Fig. 1), there exist highly localized stress con- 
centration regions near the free edges, the so-called boundary layer 
regions. This phenomenon, has been suggested to be the dominant factor 
initiating failure of some composite' laminates. 

The present thesis will analyze the boundary layer by "perturb- 
ing" the exact elasticity equations with a stretching transformation. 
Solutions to these transformed equations provide a higher order analysis 
than idealized lamination theory [1].* Hence better insight into the 
interlaminar stress behavior is obtained using the perturbation 
analysis [2]. 

1.1 REVIEW OF LITERATURE 

Bogy [3] analyzed a- bonded material containing two mutually dis- 
similar orthogonal wedges under arbitrary tractions. The stress fields 
were found to contain a mathematical singularity at the intersection of 
the interfacial plane and the loaded surface. Hein [4] studied the 
residual stresses in a two-material wedge and found similar behavior. 


*Numbers in brackets refer to- the references \isted in the 
bibl iography. 
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Hess [5] developed a plane elasticity solution for the end problem in 
a two-layer laminated strip and showed a mathematical singularity, 
defined to be a:-.point where the convergence of an eigenfunction expan- 
sion could not be 'attained, Puppo and Evensen [6] modeled the finite- 
width synmetric laminate as a set of anisotropic layers separated by 
isotropic shear layers. Each anisotropic layer was assumed to be under 
generalized plane stress, i.e., the out-of-plane normal stress 
(Fig. 1.) is zero and the in-plane stresses and displacements are the 
thickness averages of the actual values. Solutions to the corresponding 
equilibrium equations showed that while the interlaminar shear stresses 
vanish everywhere for a laminate of infinite width, they attained 
maximum finite values near the free edge of a finite width laminate. 
Furthermore, in regions far away from the free edge, the solution 
agreed well with the classical lamination theory [1]. A complete three 
dimensional analysis was carried out by Pipes and Pagano [7] using the 
finite difference technique to solve the exact elasticity equations. 

The results showed good agreements with those of Reference [6] except 
at the free edge where the interlaminar shear stress seemed to 
grow without bound for some laminates. Due to the approximate nature 
of the finite difference analysis , however, no evidence was available 
to show the intensity of the suggested singularity, Isakson and Levy 
[8] used a finite element approach to analyze a model similar to that 
of Reference [6]. Based on the constant- strain assumption within each ■ 
element, the corresponding stresses were obtained from the constitu- 
tive equations. The total-, elastic strain energy was calculated and 
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minimized [9] to yield a set of simultaneous linear algebraic equa- 
tions. Levy, et al . [10] used the same model and formulation as 
Reference -[8], -to further investigate the elastic and' plastic inter- 
laminar shear deformations in the laminate. The out-of-plane "peel 
stress" was not taken into account in both studies due to the modeling. 
Results from these solutions were quite similar to those of Reference 
[6] except at the free edge where the interlaminar shear stress txz 
was 40% -lower than that of Reference [6]. This presumably was due to 
the limitations of the finite element approximation. Improvements 
were made by Rybicki [11] who carried out a three-dimensional finite 
element analysis based on a complimentary energy formulation in terms 
of three Maxwell stress functions. These functions resulted in a set 
of simultaneous linear algebraic equations which were solved by Gauss 
reduction and the back substitution process. The "peel stress" was 
obtained in this investigation. The results showed excellent agreement 
with References [6] and [7] in regions removed from the free edges, 
while near the edges the interlaminar shear stress agreed only with 
Reference [6]; the magnitude of rxz was much lower than the singular 

, t 

value of Reference [7]. The approximate nature of the finite element 
formulation for the laminated plate apparently leads to questionable 
and quite possibly poor results at the exact free edge. Pipes [12] 
used the finite difference procedure to carry out extensive parametric 
studies including laminate geometry, fiber orientations and stacking 
sequences. The program used in Reference [12] is capable of handling 
no more than an 8 layer symmetric laminate owing to the limited 
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computer capacity. 

Several attempts have recently been made to verify experimentally 
the numerical predictions. Results by Pipes and' Daniel [13],- 
Herakovich and- Opl inger, et al...- [15] have all showed"signi.fi.cant 

stress concentration behavior near the free edges. Although stress 
intensities were not determined in these studies, there were strong 
evidences to support the numerical prediction of significant stress 
concentrations near the free edge. 

Pipes and Pagano [16] more recently developed an analytical 
solution to the elasticity equations under the assumptions of zero 
interlaminar normal stress, and zero transverse normal stress, oy, 
for the [±45]g laminate. Pagano [17] obtained yet another approximate 
solution following the cylindrical bending theory of Whitney and 
Sun [18]. Good agreement with the elasticity solution of Reference 
[7] was found for the interlaminar normal stress, (the "peel 
stress") on the midplane of a bidirectional [0/9f)]g laminate. However, 
the solution did not recognize the stress free boundary conditions 
Ty 2 (±b,z) = 0. In addition, no through thickness distribution of the 
stresses v;as available. An approximate approach was then considered 
by Tang [19] following the isotropic theory of Reiss and. Locke [20]. 

The interior domain (regions removed from the free edges) was assumed 
to be in a state of plane stress, the axial displacement u was assumed 
to be a function of x only, and the displacement components, v and w, 
were both assumed to vanish identically. The boundary layer equilib- 
rium equations coupled v;ith the compatibility equations were split into 
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two problems. Namely, a modified torsion problem and a modified 
plane strain problem. The resulting fourth order differential equa- 
tions were, solved by asymptotic expansion in terms of the ply thickness 
h/2. The matching of the boundary layer solution with the interior 
domain solution was satisfied by the imposed boundary conditions for 
the two problems. The results for a [±45]s graphite-epoxy laminate 
showed good agreement for the interior regions with References [6] and 
[7] while the interlaminar shear stress at the free edge was found 
to be lower than the predicted singularity of Reference [7]. The 
through-thickness stress distributions showed nonzero shear stresses 
Txz and on the free surfaces z = ±h as well as on the midplane 
z = 0. Also, the out-of-plane normal stress vanished on both the 
interfacial planes z = ±h/2 and the midplane z = 0. This is unlike 
the results of Reference [7] which indicated maximum values of on 
the midplane of a [0/90]g laminate and on the interfacial planes of a 
[iASjg laminate. Finally, it should be noted that the approximate 
nature of the formulation in Reference [19] did not satisfy the vanish- 
ing stress boundary conditions (±b,z) = 0 and oy (+b,z) = 0 for 
each layer. 

1.2 THE FINITE DIFFERENCE SOLUTION 

In view of the discussion in Section 1.1, the finite difference 
solution of Reference [7] seems- to serve as the most dependable solu- 
tion known to the researcher. This is due to the fact that the formu- 
lation was required to obtain the exact elasticity solution to the 
problem. However, there were inherent deficiencies in the finite 



7 


difference procedures as pointed out by Pagano and Pipes- [21]. To 
this end, numerous tests were carried out by this author to examine 
the "exactness" of the solution in Reference [7] with emphasis on Its 
behavior near the free edge. The .following observations can be' made. ’ 

(!) For bidirectional laminates {O'’ and 90" plies), all stress 
free boundary conditions were satisfied except at the four corners 
of the laminate where the out-of-plane normal stress did not vanish 
but took on a low value. Also, the sign of the outer layer 
exact free, edge was found to be inconsistent with that of the inner 
layer. These results may be attributed to the dissatisfaction of the 
equilibrium equations on the free boundaries as can be seen in the work 
of Pipes [12]. 

(2) For angle-ply laminates neither of the stress free 

boundary conditions, ay(±b,z) = T^y(±b,z) = 0, was satisfied at the 
intersection of the interface and the free edge. Both Oy and r^^y at 
this suggested singularity were of an erroneously large order of 
magnitude--as high as the axial stress ax. As a result, the inter- 
laminar shear stress tx 2 attained a maximum finite value rather than 
the possible infinity predicted by Pipes and Pagano [7]. Failure to 
satisfy the vanishing stress boundary conditions at the four corners, 
az(±b,±h) = was found again. Moreover, the sign 

reversal of stresses which was found for the bidirectional laminates 
as a result of change in the stacking sequence, was not observed for 
the angle-ply laminates. The above boundary violation may be due 
to errors inherent in the solution procedure for the* angle-ply 
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laminates. 

It may be concluded that despite the good agreement with the 
results of References [6], [11] and [19] in regions removed from the 
free edges, the .finite difference solution yields poor results near 
the mathematical singularity. In order to determine the proper order 
of magnitude of stress intensity near the singular point, a more 
rigorous analytical solution to the field equations must be obtained. 
Such a solution was described by Pagano and Pipes [21] as a "mathe- 
matical nightmare." 

The present thesis seeks a solution which predicts accurate 
interlaminar free edge stress intensities for laminates. Due to the 
above-mentioned mathematical complexities, it is certainly not an easy 
task. As described in the’ preceding section, all the previous investi- 
gations show a common result— the plane stress lamination theory is 
recovered near the central plane y = 0 provided the laminate is suf- 
ficiently wide (b/h » 1). This suggests that the boundary layer 
effect is directly related to the geometrical ratio b/h, and that the 
stress distribution throughout the laminate is the combination of the 
interior region solution and the boundary layer solution. The method 
of solution employed in the present’ thesis is the perturbation 
analyses [2, 22] developed in the 1940's to solve boundary value prob- 
lems in fluid mechanics and extended to problems in solid mechanics in 
the 1950's. The isotropic theory of Reiss and Locke [20] and the 
anisotropic theory of Tang [19] were essentially based upon such 
analyses. The main differences between the present thesis and the 
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theory of Reference [19] are summarized as follows. (1) The present 
thesis is based upon the displacement formulation in which the compati- 
bility equations are satisfied automatically. The resulting field 
equations- are' second order partial differential equations in terms of 
the displacement functions. Reference [19] was based upon the stress 
formulation in v/hich satisfaction of the compatibility equations 
resulted in fourth order^ partial differential equations in terms of 
the stress functions. (2) For the interior regions, the present 
thesis determines the three dimensional solution to the reduced govern- 
ing equations (h/b -^- 0) while satisfying the symmetry and antisymmetry 
conditions, the displacement condition, the continuity conditions 
and the vanishing stress boundary conditions on the top and bottom 
surfaces. In reference [19] the displacement components v and w 
were both assumed to vanish identically for the interior regions -and 
the axial displacement u was assumed to be a linear function of x 
alone for such regions. (3) For the boundary layer region, the 
present thesis removes mathematical complexities by considering the 
free body diagram of an infinitesimally thin slice containing the 
interfacial plane. Such a limiting analysis provides sufficiently 
accurate determination of the coefficients of the boundary layer 
solution for h/b « 1. The physical validity of the composite solution 
(interior and boundary layer solutions combined) is insured by the' 
following requirements. The material immediately adjacent to the 
Interfacial plane must satisfy the "stretched" governing differential 
equations, the matching principle of perturbation theory (Section 1.3), 
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the force and moment balance with the stress resultants on the central 
plane (y = 0), the self-equilibrating conditions of the out-of-plane- 
normal stress resultant, and the free- edge stress boundary conditions. 
Reference [19] considered two separate problems for the boundary layer 
region according to the even and odd nature of the stress components. 

A stress function following the' isotropic torsion problem was assumed 
for the T° (modified torsion) problem. A similar function was then -. 

-■V 

chosen as the. particular solution to the fourth order equation of the 
P'* (modified plane strain) problem. This particular solution along 
with the homogeneous solution (5^^ order polynomial) constituted the 
solution to this problem. The combination of the T° problem and the 
P° problem. failed to satisfy some of the stress boundary conditions at 
the free edge and on the free surfaces. And the approximate nature 
made the through-thickness stress distributions incapable of properly 
describing the force and moment equilibrium and the self equilibrating 
condition at any level of z. 

1.3 PERTURBATION METHOD 

Consider the differential equation 


-y'+y=0 , 

(1.1) 

y(0) = a , y(l) = 3 

(1.2) 


where 0 < £ « 1. 

Assuming the exact solution to the problem is not available, the 
following approximate steps must be taken: 



11 


As e vanishes, Equation (1.1) reduces to 

y' - y = 0 (1.3) 

which has a solution of the form 

y° = ae^ (1.4) 

where the superscript ° denotes the solution corresponding to c = 0 
and a is an unknown coefficient. Solution (1.4) can satisfy only one 
of the boundary conditions (1.2), For the other boundary condition 
to be satisfied, a stretching transformation is introduced in the form 

C = |B - xl/e^ (1.5) 

where X > 0 and B is the boundary limit of the stretched end (0 or 1 

in the present problem). It will be shown that this transformation 
magnifies a small region called the boundary layer In which y changes 
rapidly in order to retrieve the dropped boundary condition at the end 
X = B; Solution to the boundary layer equation must match the solution 
of the reduced equation (1.3) according to Prandtl's matching 
principle [2], ' ' 

m y ^ lim 
■ B c 

where y^^ is the boundary layer solution. •' 

It may be shown [2] for the present problem that the boundary 
layer exists near the end x = .1 and the value of x in Equation (1,5) 
is found to be 1. Hence 
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y° = Cte^ (1.7) 

since it must satisfy the first of Equations (1.2). Also, the 
stretching transformation (1.5) becomes 

5 = ^ 4 ^ ( 1 . 8 ) 

Equation (1.8) is now introduce(d to transform the original Equation 
(1,1) into 

4 + (1.9) 

for e « 1. 

Equation (1.9) has the solution 

^ c + de“^ (1.10) 

which should satisfy the second of Equation (1.2). Hence, 

c + d = 3 Jl) 

The matching principle (1.6) is now applied as 
lim y° = lim y^*- 

X -> 1 ^ ^ CO 

or, 

c = cte 

Hence from Equation (1,11)-, 

d = 3 ae (1.14) 


( 1 . 12 ) 

(1.13) 
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which yields 

yBL = + (g - ae)e"^ (1.15) 

Finally a , uniformly, valid solution is formed according ’to' the 
equation 

Yc (1.16) 

where y^. is the composite solution and (y ) represents the common 
p.art contained in both solutions. 

It is clear that in the present problem 

= lim = .lim y*" = ae , (1.17) 

C ^ “ X 1 

hence the composite solution to the original equation is 

yc = ae^ + (6 - ae)e~^ (1.18) 

The above derivation was required for the zeroth order problem 
of Equation (1.1). For a very small e, the zeroth order composite 
solution (1.18) provides sufficient accuracy. For a relatively larger 
e, solution to higher orders must be carried out to achieve better 
accuracy. This is shown in the following steps. 

. The solution to the original equation (1.1) may be expressed as 
an asymptotic expansion of the form 

GO 

y= E e\U) , e«l (1.19) 

n=0 " 

Substituting (1.19) into Equation (1,1) results in 
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2 - ^"y'n + = 0 (1*20) 


Since this is an identity equation in the nonzero parameter e, 
the coefficients .corresponding to each n must vanish for all x' in the 


domain specified by (1.1). Hence, 

y'o-yo = o 0*21) 

y‘n - yn = yn-i • " ^ 1 (1*22) 

Also, substituting (1.19) into the boundary condition (1.2) leads to 

yo(0) = a (1.23) 

yoO) = 3 (1.24) 

yn(0) = ynO) = 0 n^O ' (1.25) 


It. is cl ear that the zeroth order problem is defined by Equa- 
tions (1.21), (1.23) and (1.24). The composite solution to this 
problem can be shown to be identical to (1.18). Also, it is seen that 
at any level of approximation n, yp_-| is known, hence y^ for any n is 
■given by the first-order equation (1.22). Therefore, the stretching 
transformation (1.8) should be continually introduced near the end 
X = 1 where the boundary condition is dropped. If the asymptotic 
expansion 

CO 

y = I z\U) 

n=0 


] 


e << 1 


(1.26) 
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is assumed, the transformed equations are 


dyo ^ Q 


d? 


2 d? 


dC 


2“ - dT = -^n-1 " ^ ■■ 


And the boundary condition at x = 1 becomes 

yo(c = 0) ='s 

* 

yn(? = 0) = 0 n ^ 1 


(1.27) 


(1.28) 


(1.-29) 

(1.30) 


At this point it must be noted that Prandtl 's matching principle 
(1.6) fails to match expansions containing higher-order solutions. 
Instead, Van Dyke's matching principle [2, 22] should be employed to 
obtain a composite solution. For simplicity in the present thesis, 
only the zeroth order problem will be considered, hence no elaboration 
will be given. Nevertheless it may well be an intriguing topic. of .. 
future study. 



Chapter II 
PROBLEM FORMULATION 

F-igure 1 shows a balanced symmetric laminate of 2m plies of 

homogeneous orthotropic lamina oriented at angles [01/62/03/ /^m^s 

with the x.axis. The laminate thickness is small compared to other 
dimensions, i.e., the length dimensions are of the order L > b » h. 
One of .the orthotropic axes of the laminate coincides with the z axis. 
The laminate is subjected to a constant inplane axial strain e^. As- 
suming elastic response exists everywhere throughout the laminate, the 
field equations can be derived as indicated in the following section. 


2.1 GOVERNING FIELD EQUATIONS 

Introducing a rotational transformation {Reference [1]) to the 
layerwise orthotropic material leads to the following constitutive 
equations with respect to the reference coordinate axes xyz 



^11 ^12 ^13 ° 

^12 ^22 ^23 0 
^13 ^23 ^33 0 

000- C44 

0 0 0 C45 

Pi 6 ^26 ^3g 0 



(k) 


( 2 . 1 ) 


where the superscript k denotes the k^^ layer in the laminate . The 


strain-displacement relations in each layer are 
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•e^ = u,X 
Ey = v,y 
£2 = W,z 
Yy2 = W,y + v,z 
Yxz = W,X + U,Z 
Yxy = V,X + U,y 


( 2 . 2 ) 


where a’ comma denotes partial differentiation. 

Since the long laminate is loaded only at its ends x = ±L. 

Saint Venant's principle [23] can be invoked such that ’the stresses in 
regions far away from the ends are independent of x. Thus, the 
equilibrium equations in such regions reduce to 


3y 


^~^xz 

9Z 


= 0 


9y az 


0 


f 


3y 


+ 



(2.3) 


Combining equations .(2.1) and (2.2), and integrating the 
resulting stress-displacement relations (independent of x) results in 
the following displacement fields for each layer. 


u = (C-|y + C 2 Z + C 3 )x + U(y,z) 

V = (C 4 Z + Cs)x - Ci|-+ V(y,z) (2.4) 

2 

w = -C 4 xy .+ Cgx - Cg-- + W(y,z) 

where C-| through Cg are unknown constants and II, V and W are unknown 
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functions of y and z only. 

The following symmetry and antisymmetry conditions must be 
imposed: 

u(x,y,z) = u(x,y,-z) 
v(x,y,z) = v{x,y,-z) 

w(x,y,z) = -w(x,y,-z) (2.5) 

v(x,y,z) = -v(x,-y,z) 
w(x,y,z) = w(x,-y,z) 

V 

and the experimentally verified [13] condition 

■ u(0;y,h) = -u(0,-y,h) (2.6) 

is imposed. 

Equation (2.6) leads to a more general antisymmetry condition 

u(0,y,z) = -u(o,-y,z) (2.7) 

for ‘continuity consideration.- 

At this point, the even and odd nature of the displacements u, v, 
w in relation to y and z can readily be seen. Substitution of Equa- 
tions (2-.4) into Equations (2.5) and (2.7) results in 

'Cl = C2 = C4, = C5 = C5 = 0 (2.8) 

and 

U(y,z) = U(y,-z), V(y,z) = V(y,-z), W(y,z) = -W(y,-z) 

(2.9) 

U(y,z) - -U(-y,z), V(y,z) = -V(-y,z), W(y,z) = W(-y,z) 
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This greatly reduces the layerwise displacement field functions ( 2 . 4 ) 
to 

u = C3X + U(y,z) 

■ V = V(y,2) ( 2 J 0 ) 

w = W(y,z) 

As defined in the beginning of the present chapter, the laminate 
is subjected to a uniform axial strain. Hence the constant C3 in 
Equations ( 2 . 10 ) is nothing but the applied strain 

Combining Equations ( 2 . 1 ), ( 2 . 2 ), ( 2 , 3 ) and ( 2 . 10 ) results in 
the following set of simultaneous partial differential equations with- 
in each layer. 

C66U,yy + 05511,22- + C25V,yy + C45V,Z2 + (C35 + C45)W,yz = 0 
^ Z 6 ^»yy 0450,22 + C22V,yy + C44 V,z 2 (C23 + C44)W,yz = 0 (2,11) 

(C45 +C35)U,y2 + (C44 + C23)V,yz + C44W,yy + C33W,22 = 0 

The appropriate traction- free boundary conditions are (Fig.- I)-- 

oTy (±b,z) = 0 

Txy(±b,z) = 0 (2.12) 

Ty2(±b,z) = 0 

along the free edges, and 


cT2(y,±h) = 0 

tyz(y,±h) = 0 


( 2 . 13 ) 
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on the top and bottom surfaces of the laminate. 

Equations (2.12) and (2.13) may be expressed in terms of the 
unknown functions U, V, W in the form 

fk) 

{Ci2£x + C22V,y(±b,‘z) + C23W,z(±b,z) + C26U,y(±b,z)} = 0 

+ C26V,y(±b,z) + CggW.zftb.z) + C66U,y(±b,z)}^'^^ = 0 (2.12) 

{C44V,z(±b,z) + C44W,y(±b,z) + C45U,z(±bz)} = 0 

tCis^x + C23V,y(y,±h) + C33W,z(y,±h) + C35U,y(y,±h)}^^^ = 0 
{C44V,z(y,±h) + C44W,y(y,±h) + C45U,z(y ,+h) } ^ = 0 (2.13) 

{C45V,z(y,±h)+ C45W,y(y,±h) + C55U,z(y,±h) ^ = 0 

where the superscripts k and 1 denote the k^^ layer and the outermost 
layer (Fig. 1), respectively. Equations (2.9) also yield the following 
restrictions on the displacement fields 

{U,z(y,0)}^'^^ = 0 

{V,z(y,0)}^"’^ = 0- (2.14) 

{W(-y,0}}^'"^ = 0 

along the midplane and 

{U(0,z)}^^.^ = 0 

{V(0,z)}^'^^ = 0 (2.15) 

{W,y(0,z)}^’^^ = 0 

along the central plane and the superscript m denotes the layer 
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adjacent to the midplane (Fig. 1). 

Equations (2.11) along with Equations (2.12) - (2.15) represent 
a well defined boundary value problem. Also, it is important to note 
that these field equations were derived for individual layers. Hence 
continuity in displacements and tractions across the interfaces must 
be enforced to insure completeness of the solution. 

Equations (2.11) - (2.1.5) can be put in the dimensionless forms 


{Q66(h/b)2u,YY + Q5 bU,ZZ + Q26(b/b)2v,YY + Q45V,ZZ 
• + (Qs 6 + Q45)(^'/b)W,YZ}^'^^'= 0 

{Q26(b/b)2u.YY + Q4 bU,ZZ + Q22(h/b)2v.YY + Q44V-,ZZ 
■ + (Q 23 + Q44)(h/b)W,YZ)}^'^5 = 0 
^(Q 45 + Q36)(b/b)U,YZ + (Q44 + Q23) (h/b)V,YZ . 

+ Q 44 (h/b) 2 w,YY + Q33W,ZZ}<*^^ = 0 

{Q12£x +^»Y(±1,Z) +^,Z(±1,2) + ^,Y(±1,Z)}^^^ 
{Q16£x + ^.Y(±l ,Z) + 53^,Z(±1 ,Z) + ^J,Y(±1 ,Z)}^^^ 
{^,Z(±1,Z) +^,Y(±1,Z) + ^,Z(±1,Z)}''^^ = 


0 

0 


0 


{Qi3Ex + ^V,Y(Y,±1) + 5^,Z(Y.±1) + 5^,Y(Y,±1)}^^^ 

{5^,Z(Y,±1) +^,Y(Y,±1) +.^,Z(Y,±1)}^^) = 0 
|^,Z(Y,±1) +?^,Y(Y,±1) + ?^,Z(Y,±l)p^ = 0 


= 0 


(2.16) 


(2.1-7) 


(2.18) 
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{U,Z(Y,0)}("’) = 0 

{V,Z(Y,0)}^'^^ = 0 (2.19) 

' {W(Y,0)}("’^ = 0 

{U(0,Z)}^‘^^ = 0 

{V(0,Z)}^'^^ = 0 (2.20) 

{W,Y(0,Z)}^’^^ = 0 ■ 

where with being the largest stiffness coef- 

ij max max ^ ^ 

ficient of -the layer, Y = ^ and Z = — , the dimensionless co- • 

* D h 

ordinates, and U, V, W and their partial derivatives being dimension- 
less quantities. 


2.2 EQUILIBRIUM CONSIDERATIONS 

Before develpping the solution procedures, it will be shown that 
a close examination of the force and moment equilibrium of a section 
of the laminate will lead to significant reduction in the mathematical 
complexity. Consider the free body diagram in Fig. 2. Let 


hence,. 



k = 1 ,2,. . . ,m 





= .l . 


( 2 . 21 ) 


( 2 . 22 ) 


The force and moment equilibrium per unit length require 


rF,-0 
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NOTE : Txy , Txz (NOT SHOWN) 



FIGURE 2 FREE BODY DIAGRAM OF QUARTER YZ- PLANE 
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FIGURE 3. PARTIAL FREE BODY DIAGRAM ORQUARTER 
SECTION. 
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which gives 


dy(0,Z)hdZ = Ty 2 ^''^(Y,Z[,)bdY 


(2.23) 


and E = 0 which gives 


rl 


A 


ay(0,Z)h^(Z - Z(^)dZ = [ a2^’^^(Y.Z|^)b2YdY 

•' 0 


(2,24) 


where Zj, = 1 - E tj is the elevation of the interface in the 
j=l 

first quadrant. 

On the other hand, the force equilibrium per unit length 
requi res 

E Fx = 0 


which yields ' 


rl rl 

Txy(0,Z)hdZ = Txz(Y,Z|^)bdY 

t ' n 


(2.25) 


-At this point, an important premise must be recognized in the 
solution method of the present thesis. It has been numerically 
observed in [24] and examined in the present study that the central 
plane stresses oy^^^(0,Z) and Txy^^^(0,Z) are essentially constant 
in each layer for h/b « 1 (Figs. 2, 3). Hence it may be expressed 
mathematically that 
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2 F„ = r = 0 

y k=i ^ K 

and 

J, ^x/^>(0.Z)t|c-0 


(2.26) 


(2.27) 


for equilibrium considerations. 

Equations (2.23) through (2.27) together serve as an important 
tool to reduce the mathematical complexities in the present thesis. 
Since the material on either side of the interface Z = must 
satisfy the governing equations (2.16) and the boundary conditions 
(2.17) and (2,20), and since the interlaminar stress distributions are 
of primary interests, the boundary layer equations will be solved by 
considering only the infinitesimally thin free body diagram about this 
interface (Fig.' 4). By doing so, the boundary value problem is re- 
placed by the free body force and moment system of Fig. 4. Thus, the 
stress boundary conditions on the top and bottom surfaces, (2.18), can 
be safely ignored. This will be elaborated upon in Subsection 2.3.2. 

:2.3 PERTURBATION SOLUTION 

As described in Section 1.3, tv/o regions v/ill be considered 
separately. ‘Namely, the interior region where the solution to the 
reduced equations (e 0) satisfies boundary conditions at one end, and 
the boundary-layer region where solution to the "stretched" equations 
satisfies the boundary condition at the other end. Matching of these 
two solutions must be enforced to insure uniformity of the resulting 
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FIGURE 4. LIMITING FREE BODY DIAGRAM OFTHE-INTERFACE ■ 
Z = Zk 
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composite solution. 

• 2.3.1 THE INTERIOR REGION 

To seek a straightforaard expansion, let 

= z enUj^(k)(Y^2:) 

n=0 

.= Z e"Vn^^^{Y,Z) (2.28) 

n=0 

= Z e « 1, k = 1 ,2,3, . . . ,m 

- n=0 

where the small parameter e represents the geometrical ratio h/b. 

Substituting these expansions into Equations (2.16) and equating 
coefficients of equal powers of e to zero result in the following sets 
of equations; 

‘ {Qbb^o'.ZZ ^45Vo,Zz}^*^^ = 0 

. 1^45^0, ZZ + Q44^o,ZZ| ^ ° (2.29) 

|Q33’^^o,Zzj^‘^^ = 0 

• {Q55^^,^ + ^45^1,22 + (P36 ^45)Wo,yz}^'^^ " ^ 

{Q45^1,ZZ ^44^1,22 + (Q23 + Q44)Wo,Yz}^'^^ = 0 • 

{(Q 45 + Q36)Uq,YZ + (Q 44 + Q23)Vq,YZ + 033^1 ,zz}^*^^ = 0 


(2.30) 
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= {%6V2,YY %5^r,ZZ Q26V2,YY %V,ZZ 

+ (Q36 Q45)Vi,Yz}^'^^ " 0 

{^26^r^2 ,YY ■*■ ^45^r.ZZ ^ZZ'^r-ZJY ^44'^r,ZZ 

+ (Q 23 Q 44 )Vi,Yz}^''^ = 0 (2.31) 

{(^45 ■*■ %6)Vl,YZ (^44 + Q23 ^Vi,YZ 

+ Q44V2,YY ^ '^33Wr,Zz}^'^^ " ^ - 2 


Now the displacement conditions (2.19) and (2.20) give 


and 


(m) 

y(Y,o) = 0 



|?(Y,0) = 0 


(2.32) 

|('^)(Y,0) = 0 

n = 0,1 ,2,. . . . 


^*^^0,2) =- 0 



(^^0,Z) = 0 


(2.33) 

(^)(0,Z) = 0 

n = O', .',2, 

k = 1 ,2,. . . ,m 



■Recognizing that the boundary layer regions exist near Y = 1 1, 

1 

the stress free boundary conditions (2.17) are dropped for this 
interior region. 

The stress boundary conditions on the top and bottom surfaces, 


Equations (2.18), yield 
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{Q13"x ■^TVn,Y(Y,±U +TWn,z(Y,±l) + ^ ,y(V,±1)}^^ ^ - 0 

+-TWn,y(Y.±l) 

{TVn,z('f.il) ^^W„,Y(V.il) -^U„,z{Y,±l)}f'> = 0 

n- = 0,1 ,2, 

For the lamination theory, a = t = t = 0 

> z xz yz 

{k = 2,3,4,...,m) must also hold for the interior region. Hence 
Equations (2.34) may be generalized to 


^23 


^33 


^36 




Q13^x + — V„^y(Y,±l) + -f- Wn^zt^i') + T- Un,Y<''.il) 

{t + T Wn.Y(^i') ^ T ^n.zf''-*’)}"'’ = 0 (2.35) 

¥vn.z(''.*l) +TW„_y(Y.±l) +5|5 ^ „ 

n = 0,1 ,2, 

k = 1 ,2,. . . ,m 


The derived symmetry and antisymmetry conditions (2.9) lead to 


Un^'"^-(Y,Z) = Upf^')(Y,-Z) 
Vn^^^(Y,Z) = Vn^'"^(Y,-Z) 
Wn^'^^(Y,Z) = -W^(k)(Y,-Z) 
Un^-^')(Y,Z) = -U^('^)(-Y,Z) 
Vn^^^(Y,Z) = -V;^('')(-Y,Z) 
Wn^'^^(Y,Z) = Wn^'^)(-Y,Z) 


n = 0,1 ,2,..,. 
k = 1 ,2 m 
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Equations (2.35) may be put in the following form with respect 
to the order of e, similar to Equations (2.29) through (2,31). 


•'!44'''*''„,2{V,±1) + a45(«U„,z(V,±I) = 0 

^ Q23^"^V._i.y(Y,± 1) + ^Y,±l) 

+ Q36^^^U^-fMY.±1) = 0 

+ Q44^^^W^-i1y(Y3±1) = 0 
Q45^'^^V,Z^'^^(Y,±1) + Q55^'^^U^^2^''V,±1) 


.(2.37) 


(2.38) 


r > 1 


Thus, the interior region problem is redefined by the infinite 
sets of equations with respect to the order of e. 
the Zeroth Order Problem: 

Equations (2.29), (2.32), (2.33), (2.36) and (2-.37), (n = 0) 
consti.tute the zeroth order problem for the interior region. 

The solutions to Equations ' (2,29) have the form 


Uo^''^ = aJ*"^Y)Z + Bq^'^^Y) 

= cJ^^(Y)Z + Dq^'^^(Y) (2.39) 

= Eo(^^)(Y)Z + Fo^'')(Y) 

k = 1 ,2,.'. . ,m 


(k) 


(k) 


where ^(Y) through Fq^ ^(Y) are unknown functions. It may be 
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noted that the form of Solution (2,39) is similar to Pagano's approxi- 
mate solution of Reference Cl 7]. 

From Equations (2.36) with n = 0, it may be shown that 


The first of Equation (2.37) then leads to 
Eo'")(Y) = - 




hence 


W, 


(k) • 


^33 


W 


Z = - 


C h 

^13 ■ 


C 33 




The last two of Equations (2.37) are identically satisfied. 
From Equations (2.32) with n = 0 


Bo^’^^(O) = = 0 k = 1,2,..., 

Also, from Equations (2.36) with n = 0 


m 


B^(k)(_Y) = . Bo^*^^(Y) 
Dq^'^)'(-Y) = - Do^'^^(Y) 


{k), 


(2.40) 


(2.41) 


(2.42) 


(2.43), 


k ’= 1 ,2,. . . ,m 

-Equations (2.26) and (2.27) may now be expressed in the form 


■ k= 


FL ^23^13' 

(k)- 

L'riz ■ C 33 J 

_ 


<=22''^ Vo' 


m 

z 

k=l 


■I- E C26<'‘>h^B„''''>(Y) = 0 


(2,44) 
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m 

z 

k=l 


^16 - 


^26^13 


'33 


(k) 


h„e„b + C2j‘ 'h^D„'' '(Y) 


k x'- 


” %5<'‘)h|,Bo'<^)(Y) = 0 

+ k=1 


(2.45) 


where 8 q'^’^^{Y) and Dq'^*^^(Y) are the first derivatives of the cor- 
responding functions. Note that the -higher order terms were neglected 
in these equations. 

Equation (2.41) implies that continuity in the displacement 
W(Y,Z) can be insured only when higher order terms are included since 
(k = l,2,...,m) are different in general. 

Enforcing displacement continuity in U(Y,Z) and V(Y,Z) results 
in 

B(,^’’(Y) - B„*^>(Y) « ...... = Bo('"'(Y) '(2.46) 

= = Do'"’{V) (2.47) 

Integrating Equations (2.44) and (2.45), making use of Equations 
(2. 42) and combining the resulting Equations with Equations (2.46) and 
(2.47) lead to 


(k) 


(Y) = - 


Iq2l5 “ ^^3^3 




k = 1 ,2,. . . ;m 


(2.48) 


D 


o 


(k) 


(Y) 


( qi ^5 ' ^ 3 ^ 4 ) 
^'^2 ‘^5 “ ^ 3 ^ 3 ^ 


c,bY 


where 
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<^1 = 

% = 

As mentioned in the preceding section, higher-order approxima- 
tions are not pursued in this thesis for simplicity. Hence, the 
interior region solutions are found to-be 

u(k) ^ + 0(e) 

y(k) = + 0 (e) (2.50) 

w(k) = . ^]3 — 2 + 0 (e) k = l, 2 ,..;. ,m 
’^33^ 

where Bq^*^^(Y) and D^^*^^(Y) are given by Equations (2.48), and 0(e) 
represents the highest order term truncated in the asymptotic expan- 
sion. 

2.3.2 MODIFIED ZEROTH ORDER INTERIOR REGION SOLUTION 

Solution (2.50) does not completely satisfy the vanishing stress 
condition C2.35) to the proper order of e. This can be seen 
from Equation (2.38) where the zeroth order displacements Uq^*^^(Y,Z) 


m r 
z C 
<= 1 *- 


‘' 23^131 


12 " C 33 J 


j: 

k=l 


EC ^^^h 
E C25 \ 


(2.49) 


k=l 

m 

E 

k=l 


L^16 


^26^13' 


-33 -I 


(k) 


(k), 

.^,^66 

k=l 
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and Vq v/ere related to the undetermined first order displace- 
ment It has been described (Section 1.3) that solving 

higher order problems requires more complex mathematical procedures 
such as Van Dyke's matching principle. Also in Section 2.2, it 
ms shown that the uniform stress distribution on the central xz - 
plane (Y =0). a numerical result, is utilized as an important tool 
to reduce tnatliematical complexity for -the boundary layer region. 
Therefore, an improved zeroth order interior region solution to evalu- 
ate better stress intensity near the central plane is certainly quite 
demanding. 

To seek such an improvement, Equation (2.39) along with Equation 
(2.40) are now required to satisfy the stress conditions (2.35) 


exactly. 

Equa^: 

ions 

(2.32), (2.33) and (2,36) remain satisfied. 

The 

following 

equ5' 

ions 

are obtained. 





„ = E-O^h . 

(2.51) 





( 2 : 52 ) 




^ D,f'^Y) ' 

(2.53) 


^23 

b 

Do'(Y) ^ Bo‘(Y))^^^ = 0 

(2.54) 


where E no'- becomes an unknown constant for the layer. 

(kl 

Again enforcing continuity in displacements ^(Y,2) and 
Yo^’^^{Y,Z), respectively, yields 

= Bo^^^Y) = .... '= Bo<^)(Y) = B^[\y 


(2.55) 
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(2'.56) 


and = .... = Do^'^){Y) = Do(Y) 

Hence, Equation (2.54) may be written as 

Bo‘(Y)=0 (2.57) 


r (i^) r (i^) p (l^) 

C'3'5 Coi ‘ A 1, \ C' 




O'), 


Continuity in the displacement ^(Y,Z), as developed in Sub- 
section 2.3.1 (Equation (2.41)), will be insured only by higher order 
considerations, hence it is not imposed as a physical requirement in 
the present modification. 

Finally, recalling Equations (2.26) and (2.27) gives 


m 

E 

k=l 


^23 

Cl2 ^-F^o 


(k)^ ■ 


^k 

• 




E 

k=l 


Do'(Y) 




and 


m 

E 

k=l 


m 

E 

k=l 


r 4 - 36- p 
Cl6 + IT ^0 


Bo'(Y) = 0 


(k). 


hk + 




Do'(V) 


(2.58) 


m r ^ ^ ^ h 

T ^66 \ 

k=i ‘ » 


B ‘(Y) = 0 


(2.59) 


Since there are m + 2 simultaneous equations ((2.57) - (2.59)) 
for the m + 2 unknowns Bo'(Y), Fq'(Y), and the modified zeroth 

order interior solution can be readily determined. 

To show the improvement made in the present modified interior 
region solution, two numerical- examples are given in Tables 1 and 2. 
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TABLE 1* 



[0/90J3 (e = 
* 

0.133) 


Displacement 

Finite Difference 

ZIR** Solution 

Modified 
ZIR Solution 


0 

0 

0 

V^'^)/(£xbY) 

-0.0397 

-0.0391 

-0.0396 


-0.2467 

-0.2534 

-0.2448 

w(2)/{ej^hZ) 

-0.2055 

■ -0.2172 

-0.2072 

* Material: graphite-epoxy laminate with constant ply 

thickness. 

** Zero^b' 

order interior region. 





• TABLE 
[45/-45]g (e 

2 

- 0.133) 





Modified 

Displacement 

Finite Difference 

ZIR Solution 

ZIR Solution 

U^*<)/(ej^bY) 

0 

0 

0 

V^'^)/(ej^bY) . 

-0.7409 

-0.7298 

-0.7433 

W^l)/(Ej^hZ) 

-0,0607 

-0.2354 

-0.0604 

wf2)/(e^hZ) 

-0.0613 

-0.2354 

-0.0604 


As expected, the Modified ZIR solution yields jnore reliable 


results than the ZIR solution. 
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( 1 ) - ( 2 ) COUPLE OF Ry { Y=0) 
( 3 )- (41 COUPLE OF Rxy(Y=0) 
( 5 )- ( 6 ) COUPLE OF Rxy 


( 7 )- (8) COUPLE OF Rxy 
{9)-(IO) COUPLE OF Rxr 
(Ilh02) COUPLE OF Rz 


note; resultants of Oy AND (Tz ARE ONLY SHOWN 
IN YZ- PLANE 



FIGURE 5. FREE BODY DIAGRAM OF FIRST QUADRAMT 
OF TYPICAL SECTION 
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Hence the for,mer will be employed throughout the present thesis 
to evaluate the central plane stress intensity for the determination 

9 

of the unknown coefficients in the boundary layer region solution, 
2.3.3 THE BOUNDARY LAYER REGION 

Consider the first quadrant of the yz - plane as shown in 
Figure 5. Introducing the stretching transformation 


n 


t 1 - Y ) 

e 


(2.60) 


near the free edge Y = 1 to the governing equations (2,16) results in 
the following equations for this quarter plane of the laminate. 

QSS'J’ZZ ^ 45^ >11 ~ (^36 ^ ° ' 

^'^‘ZZ QzzV.nn Q44V>ZZ “ ^^23 + Q44^'^’nZ^^'^^ " ^ 

(k) 

{-(Q 45 + Q36)U,,z - (P 44 ^ Q2s)V3nZ + Q44W>nn + f533l'^ZZ> =0 

(2.61) 

To satisfy Prandtl 's matching principle (Section 1.3), assume the 
following composite expansions 

f [Bo(Y) + cos Og + 0(c) 

= [Do(Yj + cos + 0(e) k = l,2,...,m (2.62) 

w(l^) = [EqZ + Soe^®"' sin «o + 0(e) 

v/here B^^^^(Y), Dj^^*^^(Y) and are the Modified ZIR solution 

given by Equations (2.57) - (2.59), Po^**'\ ^re un- 

fk) 

determined coefficients, and are undetermined positive quantities 
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in radians. The subscript o implies the zeroth order solution for the 
boundary layer region. 

Substituting Equations (2.62) into Equations (2.61) and neglect- 
ing the 0(e) terms results in the following set of three simultaneous 
algebraic equations corresponding to the e° order: 

2 9^2 fkl 

^f%6^o ■ %5“o ^^0 f^26V - f^45“o )f^o “ ^^36 ^ ^45^VoSo^^ " ° 

{(Q26>o^ - Q45»o^)Po + 0)22^0^ ‘ 044“o^)'*o ' ^^23 + 1544) VoSo>‘^' = ° 

<(*'46 * ‘>36)VoPo + (<544 + t)23)Vo’’o * (544^0^ - ” 

. k = 1 ,2,. ,m (2.63) 

For each nontrivial term of Solution (2.62) to exist the 
determinants of these algebraic equations must vanish individually. 


Thus , 


%6^o^ 

- %5“o' 

^26^0^’" ^45“o^ 

■{^36 ^45^^o“o 

(k) 

^26^0^ 

" Q45“o^'. 

^22^0^ ■ '^44°fo^ 

-(Q23 + Q44)Vo 

= 0- (2:64) 

(Q45 + 

Q36?^o“o 

(Q44 + Q23)^o“o 

Q44^o^ “ ^33“o^ 



k = 1 ,2, .... ,m 


These sixth order equations may be regarded as third-order equa- 
tions. by the classical treatment [25] and the method of complex 
variables [26]. 

The six roots are found to be in the form 


{Xq( 1,2) = ± Faol^''^ 
{Xo(3,4) = ± b 
{Xo(5,6) .= ± 


(2.65) 
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where are three positive constants in terms of 

material constants of the layer (see Section 3.2). For matching 
consideration 5 however, the positive roots must 'be dropped since they 
lead to exponential growths of the displacements for large n (or small 
Y). 

Upon determining the characteristic roots from Equations (2.65), 
Solution (2.62) takes the following- general form. 

= {Bo(Y) + (P^e"®"^'^ + P2e^“o^ + 236*^“°'^) cos + 0(e) 

= {Dq(Y) + (R^e^^o’^ + + R3e^“o'i) ^ 

= {EqZ + (S-ji®“o'^ + $26^“°'' + $36^“°'^) sin + 0(e) (2.66) 

where are replaced by - 

With the above solution, the stress boundary conditions (2.17a), 
(2.17b) and (2.17c) transform to 


{[Q2g(aP-| + bP2 + CP3-) + p22(aR] + bP2 +.CR3) 

+ Q2s(S-] + $2 + S 3)]ao cos(ao Zk) 

• ■ ^23 ^22 ^26 ■ 

= - [(Qi2 ^o)^x +tDo'(±1) +-r Bo'(±l)]h} 

(2.67) 


{[Qse^^Pl ^^3^ ^26^®*^! ^^^2 ^^3) 

^36^^! ■^32 + S3)]aQ cos (cq Z^) 

- - C(Q16 + TEo)6x + ^ Oo'(t’) + 


(2.68) 
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^2 ^3) " ^2^ ■*■ S3C)] + Q45(P] + P2 + ~ 0 

k = 1 ,2,.'. . ,m (2.69) 


Note that the right hand sides of Equations (2.67) and (2.68) are all 
known quantities from the interior problem. .Since there are ten un- 
known coefficients in the k^^ layer, solving simultaneously three 
equations from the boundary conditions (2.67) through (2.69), and six 
equations from Equations (2.63) leads to the determination of the nine 
unknown coefficients in terms of gq. The accuracy of the coefficients 
thus obtained can be readily checked by the self-equilibrating condi- 
tion of the stress resultant. 


rb 


E F,.= 


= 0 


(2.70) 


for any level of Z (Fig. 2). 

Finally, imposing the moment equilibrium conditions (2.24) and 
the force equilibrium conditions (2.23) and (2.25) determines the 
values of their orders of accuracy. 

In summary, the zeroth order interior solution (ZIR) was ob- 
tained by letting h/b go to zero. The Modified ZIR solution improved 
the ZIR solution by satisfying the vanishing stress boundary condi- 
tions (2.35) exactly. The zeroth order boundary layer solution was 
obtained by transforming the governing equations and the boundary con- 
ditions (2.17) at the free edge. The matching principle was satisfied 
bythe composite solution, and the self-equilibrating condition of the 
interlaminar normal stress resultant was employed to_ check the 
accuracy of the calculated coefficients. The continuity conditions 
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in both displacements and tractions were imposed. And the force and 
moment equi'librium of the composite solution with the central plane 
stress resultants were satisfied for the interfacial plane 
Z = Zk (Figs. 3, 4). 



Chapter III 
SPECIAL LAMINATES- 

The solution method developed in the preceding chapter applies 
to balanced, symmetric laminates with variable-thickness plies. For 
certain special cases the field equations are greatly simplified by 
the vanishing of some elements in the stiffness matrix. Among the 
various laminates studied, in the literature (Section 1 . 1 ) are the bi- 
directional laminates [ 0 / 90 ]g and [ 90 / 0 ] 3 , and the angle-ply laminates 
[ 0 /“ 9 ]s C-e/ 0 ]s* These tv/o laminates will be considered in this 

chapter. 

3.1 BIDIRECTIONAL .LAMINATES WITH CONSTANT PLY THICKNESS 

When the orientation of the fibrous layer is either 0 ° or 90 “ 
with respect to the x axis, the constitutive equation reduces to 


'^x ' 

(k) 

cn 

^12 

^13 

0 

0 

0 ■ 

(k) 


ay 


Cl2 

C 22 

^23 

0 

0 

0 



02 


^13 

^23 

^33 

O' 

0 

0 


^z 

■^yz 


-0 

0 • 

0 

C44 

0 

0 



■^xz 


0 

0 

0 

0 

C 55 

0 


YXZ 

-xy- 


_ 0 

0 

0 

0 

. 0 

^ 66 . 


-^xy- 


k = 1 ,2,. . . ,m. 

Consider the laminate consisting of 2 m layers with the stacking 
sequence [ 0 / 90 / 0 / 90 / 0 / 90 . . as shown in Fig. 6a. 
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(a) BIDIRECTIONAL LAMINATE 



{b)ANGLEPLY LAMINATE 

FIGURE 6. BIDIRECTIONAL AND ANGLE -PLY 
LAMINATES 



Equations (2.16) reduce to 


^^66 

{Q22 


fhl 




fkl 

*^»YY ^»ZZ> = 0 


V»yy ■*■ Q44 ■*' ^^23 ■*' ^44) 


h 

ibj 




(3.2) 


/ ^ (k) 

{{Q44 + Q23) bl V,Y2 + Q44 [bj W,yy + Q33 ^ 


Note that the first equation is an independent equation whose complete 
solution may be assumed in the form 


u(i<) = \ 


V J(k) 

z 3p e cos otj^ Z) 
n=0 


(3.3) 


fk) 

where ' (n = 0,1,2, ) are unknown coefficients. 

Substituting Equation (3.3) into the first of Equations (3.2) 

yield 


(f<) = L Sk 
*^( 1 , 2 ) ^ /^66 h 


Hence 


|{k) = 


n = 0,1 ,2., . . . 


z a^^ sinh 


n=0 


b ^ 

1/Q66 h 


cos a,^ Z 


(k) 


(.3.4) 


(3.5) 


Solution (3.5) automatically satisfies the first of the displace- 
ment symmetry conditions (2.19) and (2.20). The second equation of the 
free edge stress boundary conditions (2.17) reduces to 
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Substituting Solution (3.5) into Equation (3,6) results in 


(k) 


hence . 


0 n .= 0,1 ,2, 

.(k) . 


U' 


= 0 ■ 


everywhere in this laminate. 


(k) 


(3.7) 

(3.8) 
(k) 


This leads to the vanishing of the shear stresses and t 

^ • xy xz 

throughout the bidirectional laminate as may be physically expected. 
The modified ZIR solution (Subsection 2.3.2) for and 
may be determined by solving Equations (2.57) - (2.59) simultaneously 
with Bo'(Y) vanishing identically. For the simplest case of the four 
layer symmetric [O/SOl^ laminate, the displacements are found to be 

(k), 


v(k) = 


D. 


o (Y) 

( 

y(2) = 


wO) = E,0)z 


(3.9) 




(lid) 


-'(2)(1) - 
(?) ^23 ^33 r 

• = FHW ° 


4- 


23 ^33 


C<Or(2)l 

^13 , ^13 ^23 

^T2T c(T)cC2) 

33 23 33 


cxh 


(3.11) 


(k) 


(Y) = 


. 0 ) 


■'-“TiTT 


i^3 " 


23 


(3.12) 


Introducing the stretching transformation 
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13) 

into- the remaining coupled equations of {3.2} results in the boundary 
layer equation in the form 

(k) 

^^22 + ^44 ^’ZZ " (^23 ^ 44 ) W»nZ> = ‘O' ^3 

{-(Q44 + Q23) V,^z + Q44 W,rin + Q33 W’ZZ^^*^^ = 0 


, (1 .yY). 
’I - h 


(3. 


Following Subsection 2.3.3, the displacements are assumed to be 




V(k) = Do^*^^(Y) + [Roe"° " cos + 0(e) 


(3.15) 


+ [Sq 6^°'’ sin + 0(e) 


v/here Do^*^^(Y), are known quantities from the modified ZIR solu- 

tion. For the four ply t0/90]s laminate, they are given by Equations 
(3.10) - (3.12). and are unknown coefficients. 

Substituting Equations (3.15) into Equations (3,14) results in 
the algebraic equations for the zeroth order boundary layer problem as 
fol 1 ows. 


,((i)22 - Q44 i>o^)Ro ■ * ‘) 44 )^o “o Sq) - 0 

{(Q44 t Q 23 )^o “0 Ro + W44 *0^ - 033 = 0 


For a' nontrivial solution, the determinant of these equations 
must vanish. Thus, 


^22 ^0^ “ ^44 “0^ "(^23 ^^44^^o "0 

(Q44 Q23^^o *^0 ^44 ^'0^ “ *^33 ‘^0^ 


(3.17) 
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(3.18) 


(3.19) 


Note that the 2x2 determinant is only a principal minor of the 
determinant in Equations (2.64). For conventional composites, such as 


graphite-epoxy and Boron-epoxy laminates. 



(3.20) 


Hence Equation (3.18) yields two pairs of real roots. For match- 
ing considerations, the positive roots are dropped, and the zeroth ’ 
order composite solution (Section 1.3) takes the following form 





(3721 y 


(3.22) 


and Prandtl's matching principle (Section 1.3) is satisfied. 
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The first and third of boundary conditions (2.17) lead to 

{[022^^1 R] + 62R2) + Q23^S]S2]cto cos (oto Z|c) 

= - £x^(Q12 + Q 22 Do'(Y) + Q23 Eo)>^''^ 

{Q44 (Rl + R2 ' SiBi - S202)>^'^^ = 0 

From Equations (3.16), additional relations between R-] 

R2 and S2> are obtained as follows 

[(Q22^1^ “ ^44)^1 “ ^^23 t ^ 

[(Q22^2^ - Q44)R2 - (Q23 Q44)62S2l^^^ == 0 

Solving Equations (3.23) 7 (3.25) simultaneously results in 

R2 ^^^j and- $2^^^ in terms/of he^/(aQ cos ao(^k - 

where 0 < c «< 1 and 2]^ is defined in Equation (2.24) and (2.25). 

Equations (2.23) and (2.24) then lead to the determination of 

and tan - ?) to their orders of accuracy. - - 

Thus the complete solution for the zeroth order displacement 

function u(k)^ are obtained and the interlaminar stresses 

between the layer and the (k+1)^^ layer, Zi, = 1 - i: ty-. , may be 

r=l 

readily calculated from the strain-displacement equations' (2.2) and the 
constitutive equations (3.1). 

3.1.1 [0/90]s GRAPHITE-EPOXY LAMINATE 

As a numerical example, the four-ply [O/OOlg graphite-epoxy lami- 
nate with constant ply thicknesses (Fig. 7a) is considered. The stiff- 
ness coefficients (after transformation) are listed below. 


(3.23) 

and S-] 

(3.24) 

(3.25) 
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(b) QUARTER Y-2 PLANE 


0<C«<l 
(c) FBD 0FZ=I/2C 




(d) [90/0]g LAMINATE 



0<C<«l 
(f)FBD 0FZ = I/2C 


FIGURE 7. FOUR PLY BIDIRECTIONAL LAMINATES 
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0 (x 10~^ psi) 90 (x 10 


c|]^ = 20.2 

c(2) 

^11 

= 0.56 

c(2) 

H2 

r(l) = 2 21 
^22 * 

p(2) 

22 

cjp = 0.56 


= 0.48 
23 

c(2) 

^23 

C0)= 2.21 

c(2) 

33 

= 0.85 

r(2) 

^44 

Cg)= 0.85 

r(2) 

■55 

o 

II 

o 

CO 

tn 



iLbiH 

2.21 

0.56 

20.2 

0.48 

0.56 

2 . 21 ' 

0.85 

0.85 

0.85 


(3.26) 


From Equation (3.8), the axial displacement function (J vanishes 
everywhere in the laminate. 

The interior region solutions (3.9) are found to be 


li 


(2) 

= Vq^ ^ = -0.0396 eyb\ 

(3.-27) 

= -0.2448 e^hZ 

(3.28) 

= -0.2072 E^hZ 

(3.29) 


From Equations (2.2), (3.1), and (3.27) - (3.29), the central 
plane (Y = 0) stresses are found to be 

ay^^^(0,Z) = 0.3552 PsO 


(3.30) 
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ay^^^{0,Z) = -0.3552 psi) (3.31) 

T^y?(0,Z) = 4y\o,Z) = o' 


Equation (3.22) gives 


= 1.10899 
02^^^ = 0.90172 
= 1.57550 
= 0.20994 


(3.32) 


Considering continuity of Equations (3.28) and (3.29) at the' 
interfaces Z = ± ~ and the exponents given by Equations (3.32), it may 
be postulated that the boundary layer effect in the 90°-ply (Layer 2) 
penetrates deeper into the interior of the laminate than that in the 
0°-ply (Layer 1). 

Hence the zeroth order composite solution (3.21) is' in the form 


= -0.0396 zyhX + [(R]- + R 2 e^^“°'^)cos ao Z] 


= -0.2448 e^hZ + [(S^ + $2 e^^“°"^)sin a'o Z] 

.= -0.0396 z^bY + [(R] + R 2 e^^“°"^)cos <xq Z] 


( 1 ) 

(1) 

( 2 ) 


(3.33) 


W. 


(2) . o Or\70 * k7 X r^C- ;6iaor, ^ , ;e2c-on)^,„ ^j(2) 


-0.2072 e^hZ + t(S-| e 


+ $26 


The unknown coefficients are found (setting Z = ;r^) to be 

• 2 
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The self-equilibrating condition j: = 0, Equation (2.70), can 

be written, in the form 

< -[ } ± C)hdn = 0 0<c.<«l (3.37) 

Substituting the coefficients of Equations (3.34) into Equations 

(3.21)' and the constitutive equation (3.1) determines the stresses on 

either side of the interface ^ be shown that Equation 

(3.37) is satisfied identically. This further confirms the correct- 
* \ 
ness of- the calculated coefficients of Equations (3.34). 

Equations (2.23) and (2.24) now become 


t)hdn = - ^3552 hcxOQ^) (3.38) 

r - T j - Ohdo = - hej^(lO^) (3.39) 

00 - 
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r - I"" he^(io^)^ (3.40) 

*^00 ^ 

1° - a^^^hn, I - dbhO - er,)drT= hexO0®)|- (3-.41) 

0 < c «< 1 

where |- is the approximate distance of the resultants 
from the interfacial plane. 

To compare with the numerical results of Pipes and Pagano [7], 
the interlaminar stresses are calculated based on the 90° ply (the 
lower layer). Equation (3.39) leads to 


^ - a (2) 


= 0,5 


0 < ? «< 1 


(3.42). 


whence = 0.180, 8.9868, 15.4505, . 

Equation (3.41) gives 

= 2.8284 


(3.43) 


(3.44) 


From Equations (3.38) - (3,41), it is clear that the stress resultant 
is of -order O(hexlO^) while the couple moment is of order 0(h^cx^ 
Hence, requiring exact satisfaction of Equation (3.39) and approximate 


satisfaction of Equation (3.41) fixes the value of ^ at 

= 8.9868 


(3.45) 


While Txy and vanish throughout the laminate, the other 
stress components are obtained in the following zeroth order form 


av^^^(n,p 0 = [2.08837 + (0.0275 e 
^ c. 


-1.8868n- 


+ 0.0051 e^^'^^®^)]£x{10^ Psi)* 


(3.46.) 
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= _ (0.3552 + 0.0546 

- 0.4098 psi) (3.47) 

\ 

“ “-3865(e’ psi) (3.48) 

?) = (0-1356 _ 0 o^gg i’ jj ,,g) 

where 0 < ? «< 1 

The last two components of stress, the 'interlaminar stresses, 
are plotted and compared with numerical results in Figs. 8-11. 

If the stacking sequence of the laminate is reversed to [OO/Olg, 
(Fig. 7d), the derivation of these interlaminar stresses is as indi- 
cated in the following subsection. 

3.1.2 [90/0]s GRAPHITE-EPOXY LAMINATE 

While .vanishes everywhere in the laminates, the modified ZIR 
solution for and can be obtained by interchanging the 
superscripts in Equations (3.27) through (3.41). 

To compare with the numerical results [7], the 0° ply is now 

* i - * 

used as the reference layer for the interlaminar stresses. 

The stress components in zeroth order forms are obtained as 


c) = [20.04 - 0.074(i®-^^®'' + JlO® psi)' (3.50) 


ay^^^(n,|-^- S) = (0.3552 + 1.5452 

- 1.9004 e®‘^°^'^)ej^(10® psi) 


(3.51) 


Ty2^^^(n',l- ?) = 7. 6996 psi) (.3.52) 


- 4) = (-1.9004 1.5452 e 

0 < C <« 1 


®-''’'*'’).,(10« psi) (3.53) 


where the last two are the interlaminar stresses. 




FIGURE 8. INTERLAMINAR SHEAR STRESS Tyz AS 

A FUNCTION OF € FOR [0/90]^. 



FIGURE 9. INTERLAMINAR NORMAL STRESS cr^. AS A FUNCTION OF'€ 
FOR [0/90]g 






FIGURE 10. PRESENT THEORY AND FINITE DIFFERENCE RESULTS FOR 

Tyz IN [O/90J3 






0 - 2 /€xX 101 psi) 



FIGURE II. PRESENT THEORY AND FINITE DIFFERENCE RESULTS 
FOR IN [0/90]s 
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FIGURE 12. 



INTERLAMINAR SHEAR STRESS Tyj- AS A 
FUNCTION OF € FOR [SO/Olg. 
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FIGURE 13. 
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FIGURE 14. PRESENT THEORY AND FINITE DIFFERENCE 






FIGURE 15. PRESENT THEORY AND FINITE DIFFERENCE RESULTS 
FOR a-z IN [90/0]g. 
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These stresses are plotted in Figs. 12 - 15. 


3.2 ANGLE PLY LAMINATES WITH CONSTANT PLY THICKNESS 

For , contemporary fiber-reinforced composites having three mutual- 
ly perpendicular planes of elastic symmetry, vanishes. If the 

laminate consists of one material with symmetric [e7-e]s of [-e/ojg 
orientations (Fig. 6b), It is called an angle ply laminate and the 
following relations between material constants are found, 

Cij^^ = , i = 1,2.3 and j = 1,2,3 

. , n = 1,2,3 * 


The modified ZIR solution gives 


Uo 

Vo 


n) = u (2) _ Q 

o • - 

0) = V ^2) = . 
0 


(C^2^33 “ ^13^23^^ Y - 
(^ 22^33 “^ 23 ^ 23 )^^^ 




^^ 13^22 ~ ^ 12 ^ 23 ^^^ 

(C 22 C 33 - C23C23)^^^ 


Z 


(3.55) 


On the central plane (Y = 0), the stresses are. obtained from 
Equations (3.55), (2.2), and (3.1) as 

c7y^^^0,Z) = - ay^2)(o,z) = 0 


(3.56) 
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Txy^^^(O.Z) = 

= [^16 


Txy^^^(0,Z) 

Cg6(C-[2^33 ~ ^13^23^ 

(^22^33 ■ <^23^23^ 


) 

(3.57) 


The first equation, indicates that the zeroth order solution 
(3.55) contributes no transverse normal stress throughout the angle -ply 
laminate. For the laminate to be in equilibrium, two self-equilibrat- 
ing conditions in addition to Equation (2.70) should be expected to 
hold (Fig. 16). Recalling Equations (3.38) through (3.41), the follow- 
ing equations may be established. 


r C)hdn = 0 k = 1,2 


0 < ? «< 1 


T ± - En,)dTi = 0 k = 1,2 

•CO 


(3.58) 


(3.59) 


Furthermore, the shear stress resultants sficl Ryw^^^ must 

r Ay Ay 

also be in equilibrium (Fig. 3 with m = 2) as indicated in the follow- 
ing equations. 

■j- Tx 2 ^^^(n,|-)hdn + = 0 where = ] Tj^^^’^^(0,Z)hdZ (3.60) 


i 


00 


= 0 where Rxy^^^ “ 


.1 


(2). 


(0,Z)hdZ (3.61) 


The characteristic equation (2.64) leads to two identical sixth 


order algebraic equations for both layers. Three positive roots to 
this equation must be dropped for matching considerations. The compos- 
ite solution will be in the form of Equation (2.66) with 0(e) 
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truncated, A numerical 'example is presented in the following sub- 
section, 

3.2.1 [45/-45]s GRAPHITE-EPOXY LAMINATE 

Consider the [45/-45J3 graphite-epoxy laminate of constant ply 

4 

thickness h/2 (Fig. 17a). The stiffness coefficients (after transfor- 
mation) are 


45(x 10"^ psi) 


= 6.745 


^1^ 


1 ) - 


5.045 


= 0.521 


0 ) 

:2 

( 1 ) 


C22 ' - 6.745 


C23‘' = 0.521 

(1) 

C33 = 2.213 

= -4.506 

16 26 


0 ) 


'36 


-0.04387 


r - r - 0 85 

M4. ■ Ss - 0*85 

Cee” = 5.33 

C = 0 

^45 ^ 


-45(x 10"^ psi) 

= 6.745 

= 5.045 

C-,3^^ = 0.521 

( 2 ) 

C22 " 6.745 

0^3^^ = 0.521 
633^^ = 2.213 

83^^^ = 0.04387 . 

544^’ ‘ =5^ - ‘>•85 

C = 5 ’33 
'"66 

= 0 
45 


The modified ZIR solution (3.55) gives 

Uo ''* = = 8 

Vq^’’^ = = -0.7433 e^bY 

= -0.0604 eyhZ 


(3.62) 





(cj QUARTER 

■ XZ- PLANE OF TOP LAYER 

A 

(d)FBD 0FZ = ^2±^ 



FIGURE 17. FOUR LAYER 



(f) QUARTER YZ- PLANE 
OF TOP LAYER 



(g) QUARTER XZ- PLANE 
OF TOP LAYER 



(h)FBD 0FZ='/2±1 


- PLY LAMINATE 
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Equation (3.57) yields x^y^^(0,Z) = = 1.154€^(10^ psi). 

The boundary layer equations (2.61) through (2.64) yield the 


algebraic equation for both layers 

- 2.5460 + 1.6337 - 0.1202 = 0 (3.63) 

which is readily transformed to 

0? - 2.5460 J + 1.6337 u - 0.1202 = 0 ’ (3.64) 

1/2 

by letting ' ' Xq = ± aQto (3.65) 

Furthermore, setting u = y - l/3(-2.5460) = y + 0.8487 (3.66) 

results in y^ - ,0.5269 y + 0.0438 = 0 (3.67) 

Let y = Pq + pQ ■ (.3.68) 


and substitute it into Equation (3.67). The resulting set of algebraic 
equations are 


Po^ ^ 

•P^Pq = 0.1756 


which give 


(3.69) 


= (0.0736) {cos (2k7T + 107.3^) + i sin(2kTr + 107.3°)} 

% (3.70) 

= (0.0736){cos(2k7T + 252.7°) + i sin(2kTr + 252.7°)} , 

k = 0,1 ,2 


By applying DeMoivre’s formula [26] and recalling Equation 
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(3.68), (3.66) and (3.65), the six roots for Equation (3.63) are 
found to be 


Xq(1 ,2) = ± 1 .2364 qq 

Xq( 3,4) = ± 0.2903 ag (3.71) 

Xo(5,6) = ± 0.9659 ag 


Hence the zeroth order composite solution (Section 1.3) takes 
the form . 


U, 


= ({p, 


* Pje 2 » 


936 ^ )C 0 S agZl 


, _ 0,7433 + {(R-ie ' ° + Rge 2 ° 


-Boa n 

+ R 36 ^ ° )C0S agZ}^ ^ 


(k) _ 


:^ 1 V ^ . -® 2 V 


Wg' ' = - 0.0604 exhZ + {{S^e ^ + $26 

. -3, an (k) 

+ $36 ^ ° )sin agZ}^ ^ 

where = 1.2364 

B 2 ^*^^ = 0.2903 
= 0.9659 


(3.72) 


Satisfying the governing equations and the boundary conditions 


leads to the following equations; 
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where 


= - 0.5871 


= 0.1707 <j) 


1 


= 1.2021 
= - 0.6309 
= - 0.1813 <t>] 


R3 


0) 


1.1897 ^1 


5-]^.^^ = 1.1358 (J).| 

$2^^^ = 0.0347 ([>] 

$3^^^ = - 1.0736 

h = 


p/^^ = 0.5871 <J>2 

= - 0.1707 (j)2 

P 3 ^^^ = - 1,'2021 ( j )2 


R 


1 


( 2 ) _ 


0.6309 (J>2 


R2^2^ = - 0.1813 <t>2 


R' 


( 2 ) - 

> 

> 

( 2 ) _ 




1.1897 (j.2 
1 . 1 358 ^2 
0.0347 ^2 


$3^2^ = - 1.0736 <\,2 


ex' 


^ ao^^^cos{ao^^hl + ?) ) 


exh 


ao^^^cos(aQ^^^(i - d ) 

0 < ? «< '1 

•■ • It can be shown that these coefficients lead to 
faction of Equations (2.70), (3.54) and (3.55). Hence 
of these coefficients is confinned. 


Equation (3.61) then leads to 

( 2 ) 


f ( 2 ) 

a 


tan 


r^“o 


a. 


w 


= 0.5 


0 < ? «< 1 


which is identical to Equation (3.42'). 


(3.73) 


(3.74) 

identical satis- 
the correctness. 


(3.75), 
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Now consider Equations (3.69). It is clear that Layer 1 
(+ 45°) and Layer 2 (- 45°) are antisyirenetric in U and symmetric in 
V and W with respect to the infinitesimal thin slice (Fig. 17d)_. 
Upon enforcing exact continuity in displacements at Z = , the 

( 2 ) 


following equation is obtained 

( 1 ) 


lim cos 


a 


0) = 


lim cos 


a 


- ot, 


(2), 


= 0 


(3.76) 


which gives 


ct, 


(D = 


(2n + 1 )ir , n = 0,1 ,2, . . . . 


a " (2n + l)ir , n = 0,1,2 


(3.77) 


Hence, 


cos 


0) 


^ (2) 

n) 

= cos 

“o , ^ (2) 

2 ^ ^ “o 

2 ‘ 5 “o 



. 


(3.78) 


for 0 < c «< 1 


where and are given in Equations (3.77). 

(k) 


( 2 ) 


Thereby 


lim 

C-K) 


tan 


i a. 


(k), 


and 


tan 


(k) 


■J 


= K 


(3.79) 


(3.80) 


0 < c «< 1 


where K is a finite large positive value. 
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At this stage, assigning any large value' for K determines the 

fa (k) 

corresponding and tan ± and hence the interlaminar 

.. ♦ . 

stresses. It may be shown that the only stresses related to 




+ c (!') 

- “o 


C are and 


. However, the latter 


vanishes identically at the free edge as required by the stress free 

boundary conditions {2,17). Hence the singular behavior is found in 
f kl 

the intersection of the free edge and the interfacial plane 
Z = . This provides a definite mathematical evidence for the pre- 

dicted singularity in Reference [7] and w.ill be further discussed in 
the following chapter. 


The interlaminar stressea re plotted in Figures 18 - 20, 

3.2.2 ”[-45/45]5 GRAPHITE-EPOXY LAMINATE 

Consider the laminate of Fig. 17e. Interchanging the super- 
scripts 1 and 2 in Equations (3.62) through (3.73) gives a composite 
solution identical to (3.72). The corresponding interlaminar stresses 
are shown in Figures 21 - 23. - 




FIGURE 18. INTERLAMINAR SHEAR STRESS Tj52 )N [45/“45] 






FIGURE 19. INTERUMINAR SHEAR STRESS Tyz IN [45/-45L 






FIGURE 20. INTERLAMINAR W0R5VJAL STRESS PV IN [45/ -45] 

^ ■ S 



figure 21. INTERLAMINAR SHEAR STRESS [ ^5/4S]g 



CO 



to 


FIGURE 22. INTERLAMINAR SHEAR STRESS Tyj IN 
[-45/45]s, 





FIGURE 23. • INTERLAMINAR NORMAL STRESS CT^ IN [-*45/^5] 


Chapter IV 

RESULTS AND DISCUSSION 


In the preceding chapter, the general method of solution of 
Chapter II was applied to the special graphite-epoxy laminates [0/90]^, 
£90/0]g, [45/-45jg and [-45/45]^, To demonstrate the capability of 
the solution the results for the,se laminates are presented and dis- 
cussed in this chapter. 

4.1 THE FOUR LAYER UNIDIRECTIONAL LAMINATES 

It has been stated in Section 1.3 that the accuracy of the pertur- 
bation solution depends upon the perturbation parameter c. That is, 
the smaller e, the better the result. This will be demonstrated in 
what follows. 

The interlaminar shear stress and the interlaminar normal 
stress 02 (the peel stress) as functions of the perturbation parameter 
G are presented, respectively, in Figures 8 and 9 for the [O/OOj^ 
laminate. From the figures, it is clear that the boundary layer width 
becomes smaller as e decreases in magnitude. (Asymptotic recovery of 
the lamination theory is implied by the incomplete domain. of ^. ) It 
should be noted that the relative extreme values of the stresses are 
finite and remain unchanged as e decreases. This indicates that the 
present theory is capable of approximating the maximum value of the 
interlaminar stress intensities for intermediate as well as small 
values of e. Also, the difference between the cases *g = 0.133 and 
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e = 0.050 is much more than that between e = 0.050 and e = 0.033. 

VJhile the curve of e = 0.033 serves as the most accurate of the three 
stress results for their corresponding e, it lends confidence to say 
that for this.[0/90]s graphite-epoxy laminate, a geometric ratio of 
0.050 (= or smaller is sufficiently small to lead to good results 
using the present method of solution. 

Numerical results obtained by this author using the finite dif- 
ference program of Pipes [12] indicate that the smallest geometric 

ratio for which that program gives physically admissible result is 
2 

e = 0-133 (~^). Below this ratio, the instability in the solution' 
does not yield satisfaction of the force equilibrium ?: Fy = 0 (Fig. 2 
and Equation (2.'26)). This may be attributed to the inherent sensi- 
tivity of the finite difference approximation to the ratio of the grid 
spacings for partial differential . equations [27]. 

Comparisons between the results of the finite difference solution 

•t • 

and the present theory are presented in Figures 10 and 11 for the case 
e = 0.133. From Figure 10 it is clear that the present theory tends to 
predict a higher maximum intensity for the interlaminar shear stress 
Ty^- The boundary layer width is approximately the same for both 
solutions. Figure 11 shows that the present theory predicts a smooth, 
continuous distribution for az which identically satisfies the self- 
equilibrating condition e = 0 {Equation (2.70)) whereas the finite ' 
difference solution yields unstable results near the free edge which 
obviously do not satisfy this equilibrium condition. In regions 
removed from the free edge, both solutions indicate asymptotic recovery 
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of the lamination theory. 

Figures 12 and 13 show the interlaminar stresses for the ,[9Q/0]s 

I 

laminate— the .reversed stacking sequence from the previous example. ■ 
From the figures, the physical validity of the present theory is con- 
firmed by the sign reversals in both and as a result of force 

and moment equilibrium (Fig. 7). Again, boundary layer shifts due to 

« 

the reduction in e are observed. The maximum stress intensities of 
Xy 2 and cT2.1n the [90/0]g laminate are found to be finite but higher 
than those in the [0/90]s laminate (Figs. 8 and 9). This is due to 
the fact that in the calculation for the [90/0]5 laminate, the 0° layer 
was employed as the reference layer. On the other hand, in the calcu- 
lation .for the [0/90]^ laminate, the 90° layer was employed as the 
reference layer. 

Comparisons between the finite difference results and the present 
theory are presented in Figures 14 and 15, The present theory again 
predicts a higher Xy 2 than the finite difference solution. Also, the 
present theory yields a more acceptable distribution for the inter- 
laminar normal stress in view of the zero stress resultant require- 
ment. . In regions removed from the free edge, the lamination theory is 
recovered asymptotically in both solutions, 

4.2 THE FOUR LAYER ANGLE-PLY LAMINATES 

Pipes and Pagano [7] pointed out that the interlaminar shear 
stress Xx 2 ' in a [45/-45]g laminate tends to grov/ without bound near 
the free edge (Section 1,1). Hence the calculated maximum intensity of 
Xj ^2 by the finite difference approximation, though higher than other 
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numerical investigations [6, 8, 11], is still very questionable. It 
was discussed in Section 1.2 that failure to satisfy some stress free 
boundary conditions were observed in the finite difference solution. 
Also, these results showed no sign reversals in Oy, and Xy^ in 
consequence of reversing the stacking sequence. 

In the present theory the mathematical evidence for the singu- 
larity in tx 2 can be shown (Subsection 3.2.1) to be in terms of 
the trignometric equation 

tan (“ 2 — " = K (4.1) 


where 0 < ? «< 1 and K is a, near-singular large number. The value of 
must satisfy equation (3.75) 



0.5 


0 < ? <« 1 


(4.2). 


Obviously, the limiting analysis of the present theory (Fig. 17d, h) 
provides no unique determination of the value of K. It is only through 

I 

experimental investigation that this value may be realistically deter- 
mined. Such an investigation should be considered as a future study. 
For the purpose of comparisons, K is taken to be 20.3713, a value that 
leads to a maximum stress Intensity within the elastic limit. 

Comparisons between the results of the finite difference solution 
and the present theory are presented in Figures 18 through 23. Figure 
18 shows the variation of the interlaminar shear stress x^^ along the 
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•) 

interface Z = ^ . The near-singular free edge intensity of the present 
theory is much higher than the finite difference result and the 
boundary layer width is much smaller. Figure 19 shows the variation 
of the interlaminar shear stress xy^. Both solutions satisfy the 
stress free boundary condition Xy^ = 0 af the free edge. The negative- 
positive variation of the present theory confirms the additional self- 
equilibrating condition 

rb 

5^ Fy = j Ty 2 dy = 0 (4.3) 

(as a result of the zeroth order vanishing of ay in the interior 
region). The finite difference solution, on the other hand, cannot 
satisfy such a condition. The erroneous of the finite dif- 
ference solution at the free edge {not shown in the figures), as 
described in Section 1.2, is believed to be caused by inherent errors. 
In Figure 20 the interlaminar normal stress 02 finite difference 

solution indicates instability near the free edge; hence, no comparison 
can be made between the two solutions in this region. Since the auto- - 

matic satisfaction of the self-equilibrating condition 

* 

rb 

E F 2 = 02 dy = 0 (4.4) 

^ 0 

has- been demonstrated by the present theory (Chapter III) and can be 
observed from the figure, and since is not proportional* to the 

*As shown in Equations (3-46) - (3.53). 
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near-singular value of K, the present theory is believed to have pre- 
dicted a more accurate maximum finite intensity of the interlaminar 
normal stress. Such a determination is most important in the delami- 
nation failure mode [21, 28] of composites. Although the moment self- 
equilibrating condition (Fig. 16) is not. directly observable from 
Figure 20, the magnitude of this couple moment can be determined as 


IM 


=/; 


hcJO^ 

dy - 0 . 0027 — — ' 


fin - lb 
length 


(4.5) 


where equals 2K, a near-singular value from Equations (4.1) and 
(4.2). Hence the self-equilibrating condition of the couple moment is 
confirmed immediately. 

When the stacki.ng sequence is reversed to [-45/45J5 (Figs. 21, 22, 
23), the interlaminar shear stress txz experiences a sign change in 

order to balance the central plane shear resultant T^^^^(0,Z)hdZ 

Z 

which also experiences a sign reversal . The sign of both ly^ and 
remain unchanged. This is in agreement with the finite difference 
results (Figs. 19 and 22, 20 and 23). For c = 0.133, the finite dif- 
ference solution predicts a small uniform cy along the central plane 
(not shown in the figures) which does not change its sign, and magnitude 
for the reversed stacking sequence. For e .<0.0133 the finite dif- 
ference solution yields erroneous results for ay due to the instability 
of the solution. The present theory exhibits no such instabilities. 

It is important to note that the interlaminar normal stress- 02 
is independent of the stacking sequence and always tensile near the 
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free edge. For both the [45/-45]g and the L-45/45J5 laminates, a 
finite maximum intensity is predicted at the exact free edge (Figs. 20, 
23). This indicates that the deiamination failure mode [21, 28] should 
always be considered for reliable design of such laminate configura- 
tions. 

It is clear that the pl^sent study has obtained improved results 
for the interlaminar behavior of the [45/-45J2 and [-45/45]5, graphite- 
epoxy laminates. Since the aforementioned self-equilibrating conditions 
were originally considered for the 2m layer angle-ply laminate (Fig. 

16), the interlaminar stress variations in any angle-ply laminate may 
be expected to be similar to those in Figures 18 through 23. 


4.3 ACCURACY AND LIMITATIONS* 

As discussed earlier, the accuracy of the present theory depends 

I* 

upon the geometric ratio e = ^ • Hence, the relative order of magni- 
tude of the individual terms in the governing equations, in relation to 
^ , should be further discussed. 

4.3.1 BIDIRECTIONAL LAMINATES 

The coupled governing differential equations for bidirectional 
laminates (Equations (3.2)) are 


{Q22(r)^'',yv + Q,4,V.2 z + (Q44 + QzsK^jW.yz = 0}"^' 

<(Q44 ^ «23)t^V,yz + Q44(f)^«.yy Q33"-ZZ = 


(4.6)' 


From these dimensionless equations, it is essential that the 
order of magnitude of the coefficients of V,yY and W,yy, V,y2 and 
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W,yz» V,zz W,zz < 0(F) < 0{1)*, respectively, in order to 
properly stretch the boundary layer region with a transformation in the 
form 


Ti = .0 ^). (4.7) 

e 

Hence, if the material properties are fixed, the geometric ratio ^ 
obviously plays the dominant role. For the graphite-epoxy laminate 
with 0.133 (Chapter III), Equations (4.6) may be transformed to 


( 0 ") 


0.046 V,YY 0.208 Wjyz ~ ^ 

0.080 V,Yz 0.007 W,yy ■*" ^ 


(90“) 


0.42 V»YY ^ >ZZ ^ 0*22 W,yz ~ 0 
0.085 V,y2 ^ 0.007 W,yy + ^^21 ~ ® 


(4.8) 

(4.9) 


It may be. observed that, for this geometric ratio, a perturba- 
tion solution using the 0° layer as the reference layer should lead to 
more accurate results. - ' ‘ 

If the' geometric ratio is now reduced to 0.050 for the same 
laminate material, Equations (4.6) become 


(0^) 


0.0065 V,yy + V,2z + 0.0782 W,yz = 0 
0.030 V,Yz + 0.001 W,yy + = 0 


( 90 =*) 


0.05 V,yy + ^jZZ ^ 0.083 W,yz ~ Q 
0.032 V,yz + 0.001 W,yy ^ ^’ZZ ” 0 


(4.10) 


(4.11) 


* e represents the approximate order of the products of Q,*.; 

and e. 
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It is clear that the order of each term relative to V, 2 z or 
W,zz shrinks as e diminishes. This means that the degree of accuracy 
of the zeroth order perturbation solution is improved by the diminish- 
ing geometric ratio. Obviously, these equations provide mathematical 
evidence to support the judicious statement made in Section 4.1— 

"A geometric ratio of 0.050 or smaller leads to sufficiently accurate 
results." 

4.3.2 ANGLE-PLY LAMINATES 

The governing differentia.l equations for the ± 45® graphite-epoxy 
laminate with ^ = 0.133 are 

0.111 d,yy + ” 0.094 V,yy - 0.007 W,yz " 0 

-0.094 U,yy + 0.140 v,yy + V,zz + 0-215 W,yz = 0 (4.12) 

-0.003 0,yz + 0.083 V,yz 0.007 Wjyy 1‘^’ZZ ~ ^ 

f 

0.111 U,yy + 0»zz 0,094 V,YY + 0.007 W,yz 0 

(-45°) i 0.094 U,yy + 0.140 V,yy + V,zz + 0.215 W,yz = 0 (4.13) 

0.003 U,Yz + 0.083 V,yz + 0.007 W,yy +. W,zz = 0 

Again, the order of magnitude of each coefficient relative to 
U,zz» V»zz and W,zz 1^^ observed. From Equations (4.12) and (4.13), 
it can be safely stated that the present theory should lead to suf- 
ficiently accurate results for 0.133 or below. Hence no further 
reduction of the geometric ratio needs be elaborated upon. 

From the above discussion, the fact that the boundary layer 
penetration becomes weaker decreases can be detected simply by 
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examining the relative orders of various terms in the governing dif- 
ferential equations. The smaller the coefficients relative to the co- 
efficient one of the weaker the boundary layer 

effect. Hence the uniform stress distribution in the central plane 
(Section 2.2) is justified by the more rapid recovery of the lamina- 
tion solution. 

It must be recalled that in the* interior region of the ’present 
theory the exact satisfaction of the vanishing stress boundary condi- 
tions on- the top and. bottom surfaces, the continuities in the inter- 
laminar stresses, and the force equilibrium in the central plane were 
inforced. Also in the interior region the exact displacement continuii 
in U and V were satisfied by the modified zeroth order interior 
region (Subsection 2.3.2). For bidirectional laminates, the slight dil 
ference found in the displacement W (Table 1) for the two layers may 
be reduced or eliminated by higher order considerations. This is 
mainly due to the differences in material properties that constitute 

the governing differential equations. For the angle-ply laminates, 

• * . * ” ■ 

the exact continuity in this displacement was found .tobe satisfied 
automatically (Table 2). 

In the boundary layer region, the bidirectional laminates again 
reveal differences in the exponential functions (Equations (3.33)) 
owing to the intractable material dissimilarities. Hence no exact 
displacement continuity in this region may be inforced for the limit- 
ing free body considered in Figure 7. Nevertheless, the satisfaction 
of the symmetry conditions, the stress boundary conditions at the free 
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edge, the force and moment equilibrium about this free body, and the 
continuity in interlaminar tractions, is believed to have attained a 
solution which is -an improvement over previously available solutions. 

On the other hand, the boundary layer region solution, for the 
angle-ply laminate exactly satisfies the symmetry conditions, the 
stress boundary conditions at the free edge, the displacement and 
stress continuity conditions, and the force and moment equilibrium re- 
quirements (Figs. 16, 17). This exactness of the present theory is 
entirely due to the favorable parametric relations 

, i = 1, 2, 3 and j = 1, 2, 3 

Ckk^^ = ’ k = 4, 5, 6 (4.14) 

■ r l . 

^n6 ~ ^n6 » ^ 1 > 2, 3 

Moreover, it is this exactness that leads to the mathematical evidence 
for the stress singularity in the interlaminar shear stress 1 ^ 2 • 

4.4 GENERALITY AND APPLICABILITY 

From the solution method developed in Chapters II and- III, it is 
clear that the detailed solution procedures of the present theory can 
be readily programmed for a computer. The simple calculation steps 
require no approximate or iterative techniques. The generality of the, 
theory can be directly applied to variable laminate configurations 
(Fig. 1) with more layers than any existing computer program can 
possibly handle. Thermal strains can be readily included through the 
constitutive equations (2.1) to determine the induced thermal stresses 



due to the temperature drop from the curing temperature of the laminate 
or for laminates in a high temperature environment. The computer 
program thus developed would be capable of predicting the interlaminar 
stress intensities between any two layers including the midplane Z = 0. 
Failure hypotheses can then be established based on the interfacial 
plane where the stress intensity reaches a relative maximum. 

Finally the generality of the solution procedures in the present 
theory can be directly applied to explore related problems such as a 
laminate with internal free edge in the form of center holes, cracks, 
etc., a laminate subject to pure bending at the ends x = ± L, time- 
dependent boundary layer effects due to cyclic loadings, and so forth. 
The important experimental determination of the material parameter K 
may also be pursued as a future research topic. 



Chapter V 
CONCLUSIONS 


In the present thesis a general method of solution for a balanced 
symmetric composite laminate subject to a uniaxial extension has been 
developed based upon a perturbation analysis of an elastic limiting 
free body containing an interfacial plane. 

In summary of the theoretical achievements of the present study, 
the following conclusions can be made. 

(1) The solution satisfies the symnetry conditions, the stress 
free boundary conditions, most continuity conditions, and 
the force- and moment equilibrium of the limiting free body. 

(2) The solution predicts smooth continupus interlaminar 
stresses with no instabilities. ‘ 

(3) The solution provides the finite maximum magnitude of the 

interlaminar normal stress laminate configura- 

. ti.ons. 

(4) For given material properties, the solution accuracy depends 

r_ 

upon the geometric ratio ^ ~ ^ • For [0/90]^ Gr/E laminate, 

e £ pleads to satisfactory results while for ,[±45]s Gr/E 
2 

laminate, e <_^ predicts satisfactory results. 

L 

(5) For all laminates with geometric ratio, 0 < ^« 1, high • 
gradient displacement, strain and stress fields are shown 
to exist near the free edge. 
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(6) The above boundary layer effect decays exponentially to 
recover the lamination solution in the interior regions. 

(7) For bidirectional laminates the axial displacement function 
U is identically zero. Hence no r^y or xxz exist in the 
laminate. 

(8) For bidirectional laminates, the interlaminar normal stress 
^2 is finite with the sign depending upon the stacking 
sequence. For example, for a [0/90]^ Gr/E laminate, a 
maximum tensile exists at the free edge v/hiTe for a’ 
■[gO/Ojg Gr/E laminate, a maximum compressive is predicted. 

(9) For angle~ply laminates, the exactness of the solution leads 
to the mathematical evidence of singular interlaminar shear 
stresses and at or near the free edge. 

(10) For angle-ply laminates, the interlaminar normal stress 02 
takes on a finite maximum tensile value at the free edge, and 
is independent of the stacking sequence. 

(11) The solution procedure can be readily programmed for a^ 
computer. Such a generalized computer program would be 
capable of predicting interlaminar stresses between any* 
two layers of a general multi-layered laminate. 

(12) The present theory suggests vaiable means for solving 
important related problems of practical interest. 
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